Generalized Christoffel-Darboux formula for classical skew-orthogonal polynomials 



1. INTRODUCTION 
A. Random matrices 

Universality of eigenvalue correlations for different random matrix ensembles and its application in real physical 
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We show that skew-orthogonal functions, defined with respect to Jacobi weight Wa,b(2;) = 
(1 — a;)"(l + x) , a, 6 > —1, including the limiting cases of Laguerre {wa{x) = x'^e~^, a > —1) 
and Gaussian weight {w{x) — e~^ ), satisfy three-term recursion relation in the quaternion space. 
From this, we derive generalized Christoffel-Darboux (GCD) formula for kernel functions arising 
in the study of the corresponding orthogonal and symplectic ensembles of random 2N x 2A'' ma- 
trices. Using the GCD formula we calculate the level-densities and prove that in the bulk of the 
- - - spectrum, under appropriate scaling, the eigenvalue correlations are universal. We also provide ev- 

[ idence to show that there exists a mapping between skew-orthogonal functions arising in the study 

, of orthogonal and symplectic ensembles of random matrices. 

^ ■ PACS numbers: 02.30.Gp, 05.45.Mt 
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I i Universality ot eigenvalue correlations tor ditterent ranaom matrix ensembles ana its application m real pnysical 
' systems has attracted both mathematicians and physicists in the last few decades [J [3, M 0, & 0, M, [m, m, E 
d . O, [l^ [l^ [l7| • From the mathematical point of view, the study of this universal behavior of eigenvalue correlations 
for various ensembles of large random matrices require (i) evaluation of certain kernel functions involving orthogonal 
(unitary ensemble), skew-orthogonal (orthogonal and symplectic ensemble) and bi-orthogonal (unitary two-matrix 
ensemble) pql yno mials ( ii) asymptotic analysis of these polynomials. The rich literature available on orthogonal 
^ ; i, [13, inl 111, III, III \m M, [i3, tH, [S [13, [HI and bi-orthogonal polynomials [H [H, H, HE [l^ corresponding to 
• different weights have contributed a lot in our understanding of unitary ensembles. 

I For orthogonal and symplectic ensembles, the problem of evaluating the kernel functions is removed by the Tracy- 
. Widom formalism, where the kernel functions are expressed in terms of that of the unitary ensemble, involving 
orthogonal polynomials [H, [H, [13, [HI [1^ ■ Results for orthogonal polynomials (OP) are then used. 

The second and perhaps a more enriching method is to evaluate the kernel functions directly in terms of skew- 
, , orthogonal polynomials/skew-orthogonal functions (SOP) and use asymptotic properties of these functions d^l- For 
this, we need to develop the theory of SOP [1 13, i, U U, M, [H, [H, HI, M, [11 so that we can have further insight into 
orthogonal and symplectic ensembles of random matrices. 
^ I In this paper, we study statistical properties of orthogonal and symplectic ensembles of random matrices with 
classical weight using SOP. To do so, we evaluate the kernel functions (which we have termed generalized Christoffel- 
Darboux sum or GCD) for the corresponding orthogonal and symplectic ensembles of random 2N x 2N matrices. As 
N oo, study of these kernel functions require a knowledge of the asymptotic behavior of the corresponding SOP. 
These can be derived by solving the Riemann Hilbert problem ^, [4^ for SOP. The other (and perhaps easier) option 
is to obtain a finite term recursion relation between SOP and OP and use known properties of the latter (also see 
[3, [3 ) to obtain asymptotic results for SOP. We use the second option. 

OP satisfy three-term recursion relation, irrespective of the weight function ^(a;) w.r.t. which they are defined. 
However recursion-relations satisfied by SOP depend on the corresponding weight. This is because these functions, 
defined in the range [xi, X2], show skew-orthonormality w.r.t. their derivatives (jl.lGp . This brings into picture terms 
like w' {x) / w{x) which gives rise to certain local behavior of the recursion-relations. For example, for SOP defined in 
the finite range [a;i,a;2], w'{x)/w{x) may have poles at the end points, which have to be dealt with in the recursion 
relations. For SOP defined in the infinite range [— cxi,cxd], this can increase the number of terms in the recursion 
relations. 

An obvious consequence of the local behavior of recursion relations is that the GCD formulae also depend on w{x). 
In this paper, we calculate the GCD sum corresponding to Jacobi weight, including the limiting cases of associated 
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Laguerre and Gaussian weight. GCD formula for weight functions with polynomial potential has been derived in 
and is mentioned in this paper for completeness. In |3|, lJ| and the authors have found compact expressions 
for SOP in terms of OP. In this paper, we use them and make further developments of the asymptotic analysis of 
these SOP. Finally, using these asymptotic results in the GCD formulae, we study the statistical properties of the 
corresponding random matrix ensembles in the bulk of the spectrum. 

We also observe certain duality property between the two families of SOP arising in the study of orthogonal and 
symplectic ensembles. (In fact, this justifies further the use of SOP to ordinary OP in studying these ensembles 
[9, 10, 33,^,^, j36, 37, 38, 39;, SO, Sll.) However, this can at best be termed in physics literature as "experimental 
observations". We do not have a clear theoretical understanding of this duality. The answer perhaps lies in the 
existence of certain ortho-symplectic group which shares such duality property. But the author's knowledge in this 
field is severely limited. 

We consider ensembles of 2N dimensional matrices H with probability distribution 

Pp,N{H)dH = -^exp[-[2TrM(i7)]]diJ, (1.1) 

where the matrix function u{H) is defined by the power expansion of the function u{z). The parameter /3 = 1, 2 
and 4 correspond to ensembles invariant under orthogonal, unitary and symplectic transformations. In this paper, we 
study P = 1 and 4 cases. The partition function is given by 



» 2N-1 

'-P.n:^ / exp[-[2Tru(i/)]]di/-(2iV)! J] 5^, (1-2) 



where is a set of all 2N x 27V real symmetric (/3 = 1) and quaternion real self dual (/3 = 4) matrices. dH is the 

standard Haar measure, g^^"^ are normalization constants for SOP [27] corresponding to P — 1 and 4. 

B. Skew-orthogonal polynomials: relevance in orthogonal and symplectic ensembles 

From the invariance of these ensembles under orthogonal (/3 — 1) and symplectic (/3 — 4) transformation, joint 
probability of eigenvalues (xi, 0:2, ... , X2n) is given by [28j : 



^ 2N 

PP,n{xi,X2,...,X2n) = \I\2n{xi,X2, ■ ■ ■ ,X2N)f W w{Xj), (1.3) 

J = l 

where A2Ar(a;i, a;2, ■ • ■ , a;2jv) = Y\j<k^^j ^ ^k) is the Vandermonde determinant. The weight function w{x) — 
exp[— u(a;)] is a non-zero and non-negative function on the interval [a, b] and have finite moments. 
The n-point correlation of eigenvalues is given by 

Rlj^\xi, . . . ,X„) = — — — / dXn+l ■■■ dx2NP0,N{xi,X2 ■ ■ ■ ,X2n), n = 1, 2, . . . . (1.4) 



(2iV-n)! 

To evaluate such integrals, the joint probability distribution Pp,N{xi, X2^ ■ ■ ■ tX2n) is written in terms of quaternion 
determinants (i.e. determinant of a matrix, each of whose element is a 2 x 2 quaternion) satisfying certain properties 
[2^. Finally, using Dyson-Mehta theorem (page 152 of [29|), one can calculate (|1.4p . For example, the two-point 

function R^f \x,y) and the level density r[^\x) is given by 

Here S is the kronecker delta. In terms of SOP (f>l(^\x) and ip^\x), to be defined in p.lOp and ()1.15p respectively, 
they are expressed as: 

2N-1 

Si^J{x,y) :^ ^ (x)^f (y)=$(«(2;)n*<''n^) 

j,k=0 2N 

= -^^^\x)Y[^<^P\y), (1.6) 

2JV 
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2N-1 

j,k=0 2N 



2N-1 

j,k=0 2N 



2N 



where 



n 

cpll'^x) ^ ^=7T(f\x)w{x), 4«(x)=^c^'')x^ /3=1,4, (1.10) 



are normalized SOP of order n. 

Here, 112^ ~ diag(l, . . . , 1, 0, . . . , 0) is a diagonal matrix and ^^^^x) and '^''^\x) are semi-infinite vectors: 



2JV 



¥'^\x)^{^i^^\x)...'^lf^^\x)...)\ ¥l^\x) = -^^f^^\x)Z, (1.12) 



with each entry a 2 x 1 matrix: 

'^i'H^) - ( ) , m^) - ( \. (1.13) 

(similar for ^''n\x) and ^\!^\x)). The anti-symmetric block-diagonal matrix Z is given by 



Z ^ (1.14) 

such that Z = -Z' and = -1. 
For 

f (-^Hx) = <i>f )(x), vl;(i)(a;) = / $(i)(y)e(a: - y)^;/, n e N, (1.15) 
these polynomials satisfy skew-orthonormal relations w.r.t. the weight function w^(x)^: 

j^^i^\xM\x)dx = 5^ra[l J), n,meN. (1.16) 
Finally, from (fTTS)) and (fL6l) . (fL7| and pTSj) we get 

^^;^(.,.)^^^^, (1.17) 

4^^(.,y) = ^%^, S^k^,y) = -'-^%^. (1.18) 



To observe the dual property among the two famiUes of polynomials 7r„ (a;), /3 = 1, 4, we skew-orthonormalize them w.r.t. w^{x) (as 
in ll.l6| l to set both the families of SOP s on an equal footing. However to study the statistical properties of symplectic ensembles only, 
this is not needed. 
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Thus a knowledge of the kernel function S2^{x,y) is enough to calculate the correlation function. In this paper, we 
will study the finite N and large N behavior of S^j^{x,y). 

Outline of the paper: 

- In section 2, we calculate GCD formulae for the kernel function S^^{x,y) (|1.6p . /? = 1 and 4, corresponding to 
different weight. 

-In section 3, we discuss the idea of duality that exists between the two families of SOP arising in the study of 
orthogonal and symplectic ensembles of random matrices. 

-In section 4, we give a brief summary of some of the relevant properties of classical OP which will be useful in our 
study of the corresponding SOP. 

- In section 5, we use results of section 2 to (i) obtain the level densities (jl.5|) for Jacobi and associated Laguerre 

orthogonal ensembles, (ii) prove that in the bulk of the spectrum, the kernel functions {x, y)/ S^j^{x, x) and hence 
the unfolded correlation functions for the above ensembles are stationary and universal. 

-In section 6, we repeat the same calculations for Jacobi symplectic and associated Laguerre symplectic ensembles. 

- Conclusion. ^ 

2. THE GENERALISED CHRISTOFFEL DARBOUX SUM 
A. Recursion Relations 

For polynomials with weight function 

w{x) = {X2 — x)"" {x ~ Xif , Xi,X2^M. (2-1) 

skew-orthogonal in the finite interval [a;i,a;2], and having finite moments, evaluation of ^^\x) and will 
involve terms like w' {x)/'w{x) which have poles at xi and X2- Hence to obtain recursion relations we expand 

\{x — Xi){x2 — x)'^^^\x)\ and \x{x — Xi){x2 — x)(^^^\x)\ in terms of SOP (^^^\x) (jl.lip and introduce semi-infinite 
matrices P^l^^ and R'^^^ such that for /3 = 4, 

{x - xi){x2 - x){^'^^\x))' = f{x){^'^^\x))' ^ P(*)$(4)(:e), (2.2) 

x{x — Xi ){X2 - a;)($(4)(a;))' = x f {x){^^^\x))' = R^^^^'^^\x). (2.3) 

For /? = 1, we get 

{x - xi){x2 ~ x)<^^^\x) = /(.T)$(i)(a;) = P^^^^^^^x), (2.4) 

x{x - xi){x2 - x)<^^^\x) = xf{x)^^^\x) = r'^^'^¥^\x). (2.5) 

Equations (12. 4p and (|2.5p are obtained by multiplying the above expansion by e(y — x) and integrating by parts. 
In this context, the Jacobi weight function is defined in the interval [—1, 1] by 

Wa,b{.x) = {l-xf{l + x)\ a>-l,6>-l, (2.6) 

where restrictions on a and b ensure that they have finite moments. 
Associated Laguerre weight function is defined in the interval [0, 00] by 

Wa{x) = x^e-^, a > -1, (2.7) 

where restriction on a ensures that they have finite moments. 
Gaussian weight function is defined in the interval [—00, 00] by 

w{x)^e-''\ (2.8) 

From here on, we will concentrate on these classical weight functions and show that the corresponding SOP satisfy 
three-term recursion relations in the 2x2 quaternion space. 



^ For random matrix ensembles with polynomial potential, the special case of Gaussian ensemble {d = 1), which is also one of the limiting 
case of Jacobi ensemble, has been worked out explicitly in Q using GCD formula, and hence not repeated in this paper. 



5 



B. Recursion relations for SOP with classical weight 

For classical weight, we expand [(/(.t)<I>(''^ (.t)] ' and [a;(/(.T)<I>('^^ (.t)] ' in terms of ^^^\x). They satisfy the following 
recursion relations: 

/(a;)«'(4)(a;) = fW$(4)^^)^ xf{x)¥*\x) = R^^^<i>'^'^Hx), (2.9) 

= xf{x)<p'^^\x) = R^^'^¥^\x). (2.10) 

where P(^) and are semi-infinite tridiagonal quaternion matrices. The semi-infinite vectors ^^^\x) and '^^^'^{x) 
are given in and (|1.12p respectively. Now, w' {x)/'w{x) and hence '^'-'^^{x) and <i>(^)(a;) has singularity at a; ± 1 

for Jacobi, and at x = for associated Laguerre. To remove them, we have 

f{x) = (l-a;2), Jacobi (2.11) 
= X, Associated Laguerre (2-12) 
= 1, Gaussian (also true for any polynomial weight). (2-13) 

In other words, SOP satisfy three-term recursion relation in the quaternion space and is given by: 

/(x)*W(x) = Pi'Ui$iti(^) + Pitl^i'H^) + Pi'Li'i'i'liW, (2.14) 

xf{x)¥^^Hx) = Ri'l^,^i%{x) + Rit^i'H^) + Ri'LMt^ii^)- (2.15) 

f{x)<^>i'Hx) = Pi^L+i<^ji(:r) + pW*«(a;)-hPi^Li*i'2i(x), (2.16) 

xf{x)^ilHx) ^ <|,+i<^ji(x)+RWvI/W(:,) + RW_^vI;W (2.17) 

where ^i!^\x) and "^i^\x) are given in p.l3p and P^-''j? and R-'^f^ are 2x2 quaternions. Equation (|2.14p - (|2.17p can 
be proved directly using the skew-orthogonal relation (|1.16p . We leave it as an exercise. In this paper, we will give an 
alternative proof by showing that the semi-infinite matrices P^^^ and R^^'^ are tridiagonal (in the quaternion sense) 
and anti-self dual. 

In terms of the elements of the quaternion matrices, ()2.14|) and (j2.15|l can be written as: 





-l(^) / 



(.r.\ ] ^ \ p(4) ,(4) ,s + p(4) p(4) 

V2n+l\-^) / \ -'T-2n+l,2n+2 ^2n+l,2ri+3 / \ f2n+?,\-^ J J \ ^2n+1.2n ^2n+l,2n+l 

I I '^2n.2n-2 ^ 

\ -'T-2n+l,2n-2 -'T'2n+l,2n-l 



(2.18) 





(2.19) 



For /? = 1, we get similar relations, where and are replaced by 5''-^-' and ^^^^ respectively. 

For the polynomial weight, the semi-infinite matrices P^'^' and i?^^-* have d quaternion bands above and below the 
diagonal 0. Thus the Gaussian {d — 1) SOP, like the Jacobi and associated Laguerre functions, satisfy three term 
recursion in the quaternion space. 

Note: Here, we would like to mention that unlike OP, the Jacobi matrix Q^'^^ coming from the relation x^^^\x) — 
Q(/3)(j)(/3) (2;)^ /3 = 1, 4, holds little importance as they do not have finite bands below the diagonal. 
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C. Proof 

To prove that the SOP corresponding to classical weight satisfy three-term recursion in the quaternion space, we 
will prove that the matrices P'^' and i?^'*^ for Jacobi weight are anti-self dual. We use the scalar products 



(2.20) 
(2.21) 



Similarly, using ^x^/'l'*' (x), i/'m^ (2^)^ and ^x^Vn*^ (a;), V'm^ (x)^ for associated Laguerre weight and ^vi'*^ (x), V'm' (a;) 
^xi/'i'^'' (a;), V'm'' (a;)^ for polynomial weight for /9 = 4 and replacing ili'^'^^ix) by (\)^^^{x) for /3 = 1, we get 



and 



where dual of a matrix is defined as 



■= -ZA'-Z. 



(2.22) 



(2.23) 



It is straightforward to see that P^^^ and R^^'' have finite bands (one in the case of SOP defined w.r.t. the classical 
weight functions) above the diagonal. Equation (|2.22p ensures that they also have the same number of bands (where 
each entry is a 2 x 2 quaternion [2^) below the diagonal. This completes the proof. 



D. Generalized Christoffel Darboux sum 



In this subsection, we generalize the results given in Q to include GCD sum for both classical weights as well as 
weight functions with polynomial potential. 

With /(y) given in (PTT|) . we use ((Ll]) and to get. 



f{y)si'^{x,y)-f{x)s!,%,x) = f{y) 



$(4)*(a:)n^n**'^(y) 



2JV 2JV 



2N 2N 



2N 2N 



2N 



¥^\x)P^^^l[^^^Hy) 



2N 



2^, 



$W(y). 



(2.24) 



Similarly, 



y/(y)'5'2^(a;,2/) - a;/(x)5^^(y,a;) = yf{y) 



^(^^\x)l[zl[¥^\y) 



2N 2N 



■xf{x) 



^^'^\x)\{z\{-i>(^\y) 



2N 2N 



2N 2N 



<f^^^\x)zzR<^^^'zz n ^ n ^'^^^ 



2N 2N 



-$(4)(a.)J|i?(4)<j,(4)(j^) 



2N 



2N 



= ¥^\x) 



(2.25) 



2ArJ 
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Combining the two, GCD formula for symplectic ensembles of random matrices with classical weight is given by 



N> 1. 



f{y){y-x) 

For the corresponding orthogonal ensembles (/3 = 1), GCD formula is derived using similar technique. From ([1 
and dmni), we get 



(2.26) 



and 



xf{x)s':l^{x,y) - yf{ii)s''^^{y,x) = xf{x) 



27V 27V 



f{y) 



¥'^\x)Y[zY[-^>W^y) 



2N 2N 



2N. 



vl/(l)(y), 



(2.27) 



-i>i'^\x)Y[zll¥'Hy) 



2N 2N 



yfiy) 



¥'^\x)l[zl[^^'Hy) 



2N 2N 



2N. 



Combining the two, the GCD formula for classical orthogonal ensembles is given by 

^('\x)[R^'\y),U2N]'^^'Hy) 



*s'27v(^' y) 



f{x){x - y) 



N >\. 



(2.28) 



(2.29) 



Here 



R^^\x) = R^'' '' - xP^P\ 13= I A. 



(2.30) 

is different for different weights. 

For example, GCD matrix for the Jacobi symplectic ensemble (including the associated Laguerre and Gaussian 
symplectic ensemble) has the following structure: 



%^^\x) 




2N. 















-R 



(4) 



'•2JV--l,2Af+l 
R2N-1,2n{x) ^^27V-2,2Ar 







N- 






i?: 



(4) 

2N-2,2N 
(4) 



(l^2Ni^)^2N-liy) - 4>2Niy)'l^2N-li^) 



-R 



(4) 







'-2Af-2.27V 
^2JV-l,2Afl'''J -''•27V-1.27V+1 








R 








(4) 

27V-l,2Ar+l 









41v-2(2;)41v+i(2^) - 41v-2(2/)41v+i(2^) 



.(4) 



-^2Ar-l,27V -^-^2^-1, 2Af 



4/tf-2(^)4w(y) - 47V-2(y)47v(^) 



(2.31) 



Similarly, for classical orthogonal ensemble, the GCD matrix has the following structure: 



27VJ 



vl;(i)(y) 



-"^2^-2, 2Ar 



'4'2Ni^)^2N-l(y) - ^2N(y)^2N~A^) 



+ r: 



(1) 

2Af-l,2Af+l 



R 



(1) 

■2Ar-l,2Af 



-yP: 



(1) 



2Af-l,27V 



V41v-2(a;)'/'27v(y) - ^2N-2iy)'4'2Nix) 



'^2N-2i^)^4N+l(y) - ^2N-2iy)'p2N+li^) 

(2.32) 



8 



3. DUALITY 



Duality between the two families of SOP arising in the study of orthogonal {/3 — I) and symplectic (/? = 4) ensembles 
of random matrices, was predicted in [Bj and In this section, we show the existence of such duality between the 
two families of SOP corresponding to classical weight, i.e. (a;) i-^- ^'(^^(x) and '^'^■^^{x) i— > $(^^(a;). For this, we 
derive recursion relations between the two families of SOP with their corresponding OP. Apart from demonstrating 
duality, this technique simplifies the derivation of asymptotic results of the SOP. 

We expand functions ^m\x) and "^m\x), m > 1, skew-orthogonal in the range [a;i,a;2], in a suitable basis of OP 
such that: 

(i) their derivatives are continuous in the range [xi,a:2] and vanish at the end-points. 

(ii) (f)m\x) and {tpm\x)) = (j>m\x) can be written as w{x)Trm\x). 

(iii) ^m\x) and \E'm^(a;) are skew-orthonormal in the range [xi,X2] w.r.t. their derivatives. 

We expand Jacobi SOP (j)m\x) and ipm\x) in terms of Jacobi OP Pj"'^^''^^^^{x), orthogonal w.r.t. the weight 
fimction W2a+i.2b+i{x) fsee l2.6p . For associated Laguerre SOP, we expand in terms of L|°^^(x), orthogonal w.r.t. the 
weight function u;2a+i(a;) (see I2.7|) while for Gaussian SOP, the basis chosen is Hj{x) orthogonal w.r.t. the weight 
. The choice of such basis ensures that conditions (i-iii) are satisfied. 



A. Jacobi SOP 



Jacobi SOP corresponding to orthogonal and symplectic ensembles are given below. We see that there exists a 
relation between {x) and (j)"m {x) and their derivatives. We will use the following identities: 

(1 - X^)Wa^b{x) = Wa+l,b+l{x), Wa,b{.x)Wa+l,b+l{x) = W2a+1.26+1 (a:), (3.1) 

For 13 = 1, with (^'^ {x) and "^'^{x) satisfying conditions (i-iii), we have for m > 1, 

(52m) i^27yi+i{^) = Wa+i,b+i(.x)P^„^ ' (a;), (3.2) 



(52m) 02m+l(a^) = Wa,b(x)[A2m+l^'2m+i (x) - B2m-l^'2m-l' (3-3) 

(5ri)''V^2(-) - ^^^^^^^P^Pr^tY'^\-) + I'^-lJ^-.i-l (3.4) 

^2m 

(5^2)'^'4'i(^) = w^,,{x)Pl':;^'^''+'{x), (3.5) 



I52m. 

with (,g^ ') % {x) = J e{x - y)wa^{y)dy and 



52^ = 52m+l = ^2m^^'^'''^\ m = 0,l,.... (3.6) 



■1+ 

For /3 = 4, with $m-'(x) and "^\m'{x) satisfying conditions (i-iii), we have for m > 1, 

(52m) ' 4m+l(2;) = Wa+l^b+l{x)P^^}{^''"'^{x), (3.7) 
{gtli'^i^'il+M = Wa.b{x)[A2mP'2:^'^^'^\x) - B.^^.P^^+Y'+^x)], (3.8) 

i9it)''^'^li^) - ~^^^^^^PttY'^\x) + 7^^r/Vri^2(x), (3.9) 

^2m-l 

(5^t^)'^'4tl(2:) = ~Wa.bi^)P^:.tY'^\^), (3.10) 



with 



52m — y2m+l — 'hm-1 ' m — 1, Z, . . . . l^.-L-Lj 



where 



0) ^ ^ (j+2a + 2)(j + 26 + 2) j(j+2a + 26 + 2) (j + 2a + 2)(j + 26 + 2) 

^1-^-^' (j-H2)(j + 2a + 2& + 4)' ' (2j + 2a + 26+l)' ' (2j + 2a + 26 + 5) ' 
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2a+1.2fc+l ^ , 2a+1.2fc+l^ = 0, 1=1,2,.... (3.13) 

(/)'^\x) and can be calculated using (|6.ip and Gram-Schmidt method for SOP. Here we note that SOP for 

/3 = 4 is lower than that of /3 = 1 by an order 1. 

B. Associated Laguerre SOP 

Associated Laguerre ensembles of random matrices can and does play a significant role in describing real physical 
systems There exists a simple duality relation between the SOP "^m\x) and ^m\x), their derivatives and the 
normalization constant for to > 1: 

M'Jix) = -as<i>ifix); <i>l^Jix) = ~a^¥^\xy, g^Ji^ = gl\\ where '^s - ( J ) (3.14) 
is the Pauli matrix.'^ 

We now present the results for the SOP corresponding to associated Laguerre weight. We use 

2'^°'+^Wa{x)Wa+l{x) = W2a+l(2x), Wa+l{x) = XWaix), (3.15) 

the latter vanishing at x = for all a > —1. 

1/2 /-i^ /i-i 

For (3=1, with (gin ) ^m {x) = Wa{x)iim {x), we have for to > 1, 

= r+"^Wa+i{x)Lllt\2x) (3.16) 

(.9^:^+i)'^'4lUi(^) = 2"+i/^z.a(x)[AL+ii^:;+\(2^) - BL-iLltt^C^x)], B\ = 0, (3.17) 

i9i'y4'l(-) - r+y''-^^^L'^rt\{x) + 7lr-^24:^-2(^), rn^O, (3.18) 

5^™ = 32™+i = /i2™ \ m = 0,l,.... (3.20) 
where, {gi'^f'4'\x) = 2«+V2 ^(x - 2/)i/;a(y)d2/. 

1 /2 

For /3 = 4, with (g^-*) ^ln{x) = Wa{x)T:m\x), we have 

= 2'^+3/2z«,+i(x)i2^+i(2x) (3.21) 

(ff^tUi)'^'4tl+i(^) = 2'^+V2^,(x)[AL+ii^:5i\(2^) - S2^™-ii2™t\(2a;)] (3.22) 



(3.23) 



iai'jf'Mtix) - -2'^+3/2^W^^2a+l^(2^) ^ 7fc^Vri-2(-), - ^ 0, 

^2m 

= -2^+'^'wa{x)Ll'}^\2x) (3.24) 
.4tl-.9^tUi=^2::«+\ m = 0,l,.... (3.25) 



where (^q'^'') ^ 0Q'*''(a;) oc Wa{x). 

Here 7j'i^2 =7/"' ^j" ^'^'^ given by 



ij- = ^^^^ A] = = 1 + 2« + 2, (3.26) 



-^J-l , 2a+l 



^ Using 113.141 1 in 111.21 1 we can see that the partition functions corresponding to Laguerre weight also share a duality relation. 
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Here, the identity 

^{wa+ii^)Lf''+'\2x)} = ^w4x){Af^,Lf,\+'\2x) - B^_,Lf_^+'\2x)}, (3.28) 
is used to derive (|3.16p and p.lSp from (|3.17p and p.l9p respectively. 



Gaussian SOP 



For Gaussian weight, with w{x) = e ^ and /3 = 1, there exists a duahty relation between the SOP (x) and 
^rn{x) and their derivatives for m > 1, and is given by 

(9i]l+i)'^'4ll+ii^) = H^)H2m{x) (3.29) 

(322+i)'^'4™+i(^) = ^«(.T)h(l/2)i/2„,+i(x) + 2mF2™-i(a:)] (3.30) 

(ff22)'^V2™(-^) = -27«(a;)i72„,_i(a;) + 2(2™ - l)^^^l_^{x), m ^ 0. (3.31) 

(32^) '^'4™ (^) = ^(a;)i?2™(a;) (3.32) 

52m = 52m+l = /i2m, m = 0, 1,.... (3.33) 

1 /2 

with {g^^) ipQ^\x) — e(a; — y)w(y)dy and Hj{x) the ordinary Hermite polynomials. 
For /3 = 4, we have 

1/2 
) ^ 

= 2u;(x)il2™(a:) + 4m4ti-2(^) (3-36) 

(5^m)'^'4'i(^) = -u;(x)il2™+i(a;) (3.37) 

52m = 32m+l = 'iam+l, TO = 0,1,.... (3.38) 

Thus for the Gaussian weight, SOP for /3 = 4 is higher than that of /3 = 1 by an order 1. This is exactly the opposite 
of what we saw for Jacobi SOP. 



(5^m+l) 4t+li^)=Mx)H2rr,+ lix) (3.34) 
(5^m+l)'^V^tl+l(a;) = w{x)[-{l/2)H2r.^+2{x) + (2to + l)H2m{x)] (3.35) 



D. Proof 

The Jacobi SOP ^j^^ {x) and ■tp^^'' (x) can be written as 

(5fV'%r(-) = 7ri^^-a+i,+iPf_t^^^^+^(x) +x:7ri^^(5rv'%r^(-), ^ > 2, 7ir,^) ^ 0, (3.39) 

fe=0 

= 7f-.+i..+i^-r'"+^(-)+E7l^'(5i^Wl^^(.), J>1, 7?V0, (3.40) 

fc=0 

such that they satisfy conditions (i-iii). For example, Wa+i.t+iix) is used instead oiwa,b{x) in order to satisfy condition 
(i). The order of the OP are fixed by (ii), while (iii) fixes 7^"'''. 
Differentiating and using the identity 

d r f \ j-,2a+l,2h+l / \i t \i A r,2a+l,2b+l/ \ j-, r,2a+l,2b+l/ \i foAT\ 

— {Wa+l,b+l{x)P^ {X}} = Wa,b{x){Aj+iPj^^ • {x) - Bj_iP^_^ • {x}} , (3.41) 
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we get 

fc=0 



1/2 

fc=0 

1/2 

fc=0 

In this paper we will give the proof for [3 — A. The proof for [3 — 1 follows the same line of logic and can be found 
in and % 

Using p.ip and orthonormality of P^"'^^''^^^^ [x) w.r.t. the weight function 1^20+1, 2b+i (a;), the scalar products 

(4'i(^),^^'i-2.(^)) =0,^7^-1-1^=0, fc = l,2,.... (3.44) 

We also have 

(4tl(^),^^2-2fe+i(^)) = 0,=^7^r_-2lU=0, fc = 2,3,..(3,.45) 
(4tUi(^),4tL2fe(^)) = 0,=^7f:"^2.=0, A: = 1,2,..., (3.46) 
(4'i+i(^),4™-2fe+i(^)) =0,=^7^™-2.-i=0, fc = l,2,.... (3.47) 

Since odd SOP is arbitrary to the addition of any multiple of the lower even SOP, we set 72m'-i = 0. Choosing 

(2m-l) 1 , (2m) -. , 

72m-i = -A:i:zr^ and 7^„ ^ = 1, we get 

52m = 52m+l = ^2m^-i^'''*'^ • (3.48) 

Here we have used p.l3p . Finally, 

(4tl(-), 4:Li(-)) = 0, =^ 7^-3^^ = - > 2. (3.49) 

V / /i9m — 1 



Thus we get J3. 7113. lOp . 

To prove (|3.2HI3.24|) . we start with the expansion 



{9f)'"c^f{-) - 7jiV'2'^+^/^-a+i(x)X-|^(2x) +g7ri^^(5f )'''4'H-). J > 1> 7]r/^ ^ 0,(3.50) 

fc=0 

(5f V^Vf ^(a:) = 7?'^2'^+3/^«;.+i(x)L^2:^t^(2a^) +E7i^'^(5i'V^Vi^^(x), J > 1, 7^ 0. (3.51) 

k=0 

Using p.28p . we get ip^^^^x) and (f>^j'\x). Finally, we follow the same procedure and use (|3.26p - (|3.27p to obtain 
For Gaussian SOP, we expand 

i9f)''^fi^) = 7jiV^-(-)i^.(-)+E7ri^^(3f)'''4'^(-), J=0,1,..., 7ii-/^^0, (3.52) 

fc=0 

j-i 

(5fV^Vf^(x) ^ j^'^w{x)H,^^{x)+Y,ji'\gi'^)'^'i,i'\x), j>l, 7?V0, (3.53) 

fe=0 

and use the relation 



(3.54) 
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4. CLASSICAL ORTHOGONAL POLYNOMIALS AND SOME RELEVANT FORMULA 



Orthogonal polynomials Pj (x) of order j, associated with weight function w{x) in the interval [xi , X2] is defined as 



/•X2 J 

/ Pj{x)Pk{x)w{x)dx = hjSj,k, P]{x) = fc^ 

"'^1 1=0 



(4.1) 



(7) . 

where hj is the normalization constant and fc is the leading coefficient. They satisfy three-term recursion relation 



xPj{x) = Qj,j+iPj+i{x) + QjjPjix) + Qj.j-iPj-i{x), j = 0, 1, 

where Qj,k is the recursion coefficient. 

For classical weight functions (|2.6l 12.71 12. 8p , defined in the interval [xi , X2] , we have 



(4.2) 



Table 1 



Jacobi weight: [—1,1] 



Associated Laguerre weight: [0,oo] 



Gaussian weight: [— cx3,cx3] 



Pojx) 



^{a + b + 2)x + ^{a-b) 



-x + a + 1 



2x 



2"+''+^ r(j+a+i)r(j+b+i) 

(2j+a+b+l) r(i+l)r(j+a+b+l) 



a _ rU+a+1) 



, a,fc _ j_ / 2j+a+b 



2j(j+a+b) 
(2j+a+b)(2j+a+b-l)' 



_ Ml! 



2^ 



Qj- 



J-1 j' 



1/2 



(2j+a+b)(23+a+b-2) 



2(j+a-l)(j+b-l) 
(2j+a+b-l)(2j+a+b-2) 



2j + a - 1 



-J + a - 1 



J-1 



A. Asymptotic formula 



Here we give a brief summary of the asymptotic results of OP with classical weight which will be useful in our 
analysis of the corresponding SOP. 

Jacobi polynomial pj°'''''(a;), for large j, arbitrary real a, b and fixed positive number e, is written as 

x = cos 0, £ < d < TT — e, 



, a,bx-l/2 



-— - COS + ^ ^ - a + d + 0{j-'), 

TTsmfc' 2 2 2 



(4.3) 



where we have used [l^ and h'^'^ ~ 2'^'^+'')j'~^. 

Associated Laguerre polynomial Lj{x), for large j, has formula of the Plancherel-Rotach type [l^ given below. For 
'a' arbitrary and real, e a fixed positive number, we have for 



x = (4j + 2a + 2) cos^ 6*, e < 6* < 7r/2 - ej" ^ ^ 



(-!)■'" 



•\/27rj'sin0cos^? \ x 



a/2 



Svr. 



sin[(j + (a + l)/2)(sin20 -26) + —]} + {jxy^O{l), (4.4) 



where hj ~ Unlike the Jacobi case, for given x, 9 depends on j; for example 9j — 0jzpi ~ ±(2j tan^j) ^. Then 



with 9 = 6j, we can also write 



(-1) 



J-1 



-\/27rj sin 9 cos 6* \x 



. ^ a/2 



sin[(j + (a + l)/2)(sin26l - 261) ± 29Aj + —]} + {jx) ^O(l). 



(4.5) 
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Equation (|4.5p will be useful in deriving the asymptotic formulaefor SOP. 

Finally, Hermite polynomial Hj{x), for large j, has formula of the Plancherel-Rotach type given below. For e a 
fixed positive number, we have 

X ^ {2j + COS0, e<0<TT-e, 

e-^'/'H,{x) = ^jdL (-) "\h{j/2 + l/4)(sin2^ - 20) + ^] + 0{r'). (4.6) 

Unlike the Jacobi case, for given x, depends on j; for example 0j — 0j^i — ±(2jtan6'j) ^. Then with = 0j, we 
can also write 

e-^'/'H,^^ix) = ^^'^''1 (-) sin[(j/2 + 1/4) (sin 20 -20)±0+^-^]+ 0{r'). (4.7) 
where we have used again 0j — 0ji^i — ±(2j tan^j)^^. 

5. UNIVERSALITY FOR ORTHOGONAL ENSEMBLES 

In this section, we use GCD formula for orthogonal ensembles of random matrices ()2.32|) along with the asymptotic 

results of SOP [3, 0] to prove that with proper scaling, the kernel function {s'^^{x,y) / s'^^{x,x)) and hence the 
correlation function for Jacobi orthogonal and associated Laguerre orthogonal ensemble is stationary and universal. 

A. Jacobi SOP 

The asymptotic results for the SOP are derived using (|4.3p in (|3.2j|3.5p . Here, Aj « Bj « — j/2 and 7j « 1 for large 
j. For 

X — cos 0, € < < n — e, 



= \l ^ cos[(2m + a + b+^)0-{2a+^)^] + 0{2m)-\ (5.1) 



,,(1)/ ' 1 



sm( 



sin[(2m + a + b +^)0 - {2a + ^)|] + 0(1) + 0(2m) ^ 



2 

4m+i(2:) = 2m 



TTSmf 



sin[(2m + a + 6+ ^)6I- (2a+ ^)^] +0(2m) ^ 



(5.2) 



(5.3) 



41^(2^) - \/^37T cos[(2m + a + h+\)0-{2a+ |) + 0{2m)-\ (5.4) 
TT sm 2 2 2 



Equations ([ETjl . and are obtained by using (1131) in ([331) and ([XS)) . Equation is obtained by 

partial integration of (|5.4p . Here, we note that by directly differentiating ip^^^x), we can get (p'^^^x) to the leading 
order, thereby confirming our result. 

B. Universality in Jacobi orthogonal ensemble 

In this subsection, we calculate the level density of Jacobi orthogonal ensemble and show that in the bulk of the 
spectrum, the scaled or "unfolded" [11] kernel function {S'^^{x,y) / s'"2}]{x,x)) is stationary and universal. 
To study the kernel function (|2.32p . we calculate ^"2]^-! 2A" ^2N-i 2N+v ^2N~i 2W ^2N-2 2N- expand 



2m+2 
j=2m-2 
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We get 



(1) _ 



2m+3 

E 

j=2m-2 



2m+2 

E 

2m-2 



(5.6) 



(1) 



2m+l,2m+2 



\ 



52^+2 (2to + 1)(2to + 2) (2m + 2a + 26 + 3) (2m + 2a + 26 + 4) 
(4m + 2a + 26 + 3)(4m + 2a + 26 + 5) ' 



(5.7) 



d(1) _ o 



52m+2 (2™ + l)(2m + 2) (2m + 2a + 26 + 3) (2m + 2a + 26 + 4) 
M 52m + 2a + 26 + 3)(4m + 2a + 26 + 5)(4m + 2a + 26 + 6) ' 



(5.8) 



(1) 

2m+l,2m+2 



0, 



52m+2 [(26+1)^ - (2a+ l)^](2m + l)(2m + 2)(2m, + 2a + 26 + 3)(2m + 2a + 26 + 4) 



52 



(1) (4m + 2a + 26 + 3)(4m + 2a + 26 + 4)(4m + 2a + 26 + 5)(4m + 2a + 26 + 6) ' 



for 



a = 6, 



(5.9) 



R 



(1) 

2m,2m+2 



= -2. 



52^+2 (2to + l)(2m + 2) (2m + 2a + 26 + 3) (2m + 2a + 26 + 4) 
M (4m + 2a + 26 + 3)(4to + 2a + 26 + 4) (4m + 2a + 26 + 5) ' 



(5.10) 



For m = N — 1, large N, we have (52^/52/^-2) — ^ ^'^^ 



^2Ll,2iV~^'+0(7V)> 



(1) 



N 



0(1), 



^2Ll,27V^(&-«)0(l), 



i?: 



(1) 



N 



0(1). 



■•2^-1,2^+1 2 -"-2^-1, 27V V" "-;>-^V-^;, -"-2Af-2,2W 2 

Finally using (|5.1[) and (|5.2|) . defined in the interval [e, tt — e], and (|5.11|) in the GCD formula (|2.32p . we get for 

x^cosS, y ^ X + Ax = cos{9 + A9), x — y~A9sm6, 



(5.11) 



ix-y){l-X^)Si]l{x,y) = R2N-2,2N ^^iNix)^2N~liy) - ^2Niy)^2N-li^ 



+R 



(1) 

'2N-1,2N+1 
''2N-1,2N 



nN-2ixmN+iiy) - w2N-2iymN+i(^) 

^^2N-2(^)^2Niy) - ^2N~2iy)^2Ni^) 



3(1) 

2Af-l,2Af 



— [sin(/2Ar(0)) COs(/2Ar_2(^? + AO)) - sin(/2Ar(^? + AO)) C0s(/2W-2(e))] 

Zn 

+ 7^ [Mf2N-2iO)) COsif2N{0 + AO)) - sin(/2Ar_2 (^? + AO)) C0s(/2Ar ((?))] 
ZTT 



COS0 



sin(/2w-2(e)) sin(/2w(^ + A0)) - sin(/2w_2(0 + Ad)) sin(/2jv(^))] , 



2tt sin 

where, in the second step, we have dropped 0{1/N'^) term. Thus we get 

{x -y){l~ x^)S^^^{x,y) = [sin[/2Ar(0) - {f2N-2{0 + AO)] - s\n[hN-2{0) - f2N{0 + A6)]] 



2tt 



cose'sin(26') 
- 2.sin^ 



1 

'2^ 

cos sin(2 



sin[A6i/^^(6i) - 20] + sm[A9f^j^{e) + 20] 
sin[A0/^^( 



2tt sin 



■ sin[2iVA6i] [cos(26') - cos^ 



(5.12) 
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which gives us finally 



^2N (.•''■: y) 



sm{2NA9) 
7rA6lsin6i ' 

sin(2iV(l - x^y^^^Ax) 
ttAx 



With a: — > y, we get the level density p{x) 



S2Nix,x) := p{x) 



2N 



\x\ < 1. 



\x\ < 1. 



With Ax — > 0, i.e. in the bulk of the spectrum, we get the "universal" sine-kernel 

r = p{x)Ax. 



(5.13) 
(5.14) 

(5.15) 
(5.16) 



C. Associated Laguerre SOP 

To obtain the asymptotic properties of associated Laguerre SOP, we use the relations 2'^^^^'^y/xWa{x) = ■\/?i'2a+i(2a;) 
and 2°-^^^^Wa+i{x) — Vx-\/w2a+i(2x). Replacing this in (|3.16p - (|3.19p and using (|4.4p we get the asymptotic formula. 
For arbitrary a and e a fixed positive number, 



with 



2x ^ (8m + 4a + 4) cos' 
1 



AmVTT sin COS' 
2 



sin(/2™((?)) + 



0(1) 



72m, 



V2^ 



V TT tan 9 



Mf2mm + 



\/2mx 



sm f cos ( 



COs(/2™(0))+O(l) 



\j2mx 



4:^.rw = 2/^^^' 



4m±2 (2^) 



COs(/2m(^')) 



\/2m; 



1 



sin(/2™(0)T40) 



0(1) 



4m-\/7r sinflcos^ 

Similar relations hold for ((57T5)) . (|5?Tgi) and (jOO)) as m ^ m ± 1. Here 

/2™(6') = (2m + a+l)(sin26'-2e') + ^. 



V2^ 



(5.17) 
(5.18) 

(5.19) 
(5.20) 

(5.21) 
(5.22) 



D. Universality in associated Laguerre orthogonal ensemble 

In this subsection, we calculate the level density of associated Laguerre orthogonal ensemble and show that in the 
bulk of the spectrum, the scaled or "unfolded" [2^ kernel function (s'^^^ix^y) j s')2f^(x^x)) is stationary and universal. 

To calculate the kernel function ()2.32p . we need i?27V-i 2Af+iJ -^27V-i 27Vi ^27V-2 2Af ^^'^ 2N- F'or this, we 
expand 

2m+2 

= E P2ll+iA^i^l (5-23) 

j=2m-2 

2m+2 2m+3 

^ E 4'„l,V'f^(^), x'41+i{x)= E Rill+iA'^ix), (5.24) 

j=2m-2 j=2m-2 
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Using (|3T7l) in ((5?23l) . we get 



p(i) _ 

-^2m+l,2m+2 ^ 2 



y2m+2 
y2m 



(2r7i+ 1)(2to + 2). 



(5.25) 



Using ([XTO)) and ((HTT)) in (jOi)) . we get 



^2m,2)n+2 — -'i-2m+l,2m+3 



\ 



.'/2m+2 



(m + l)(2m+ 1), 



(5.26) 



and 



(1) 

2m+l,2m+2 



y2m+2 



92 



(1) 



(to + l)(2m + l)(4m + 2a + 4). 



(5.27) 



For TO = — 1, large N, we have (.92Jv/527V-2) — ^ 



P. 



(1) 



2N-1.2N 



2N^ + 0{N), 



R 



•2Af- 



2Ar 



-A^2 + 0(iV). (5.28) 



Finally using ((ET5)l . (|ETg)l and ([OT|) . defined in the 6 interval [e,7r/2 - eTV-i/^], and ([g?^ in the GCD formula 
(12:32)) . we get for 

X = (4A^ + 2a + 2) cos^ 9, y = x + Ax, 
xix-y)S^^]^{x,y) = ^2W-2,2Ar V4N(a;)V'2j^~i(j^) - V'2j^(y)V'2W-i(2^) 



-^27V-l,2Ar+l 
\^2N-1,2N 

N 



V4jV-2(2;)V'2l^+l(j^) - V'2lv-2(2/)V'2!^+l(2;) 
xP2N-U2n) ^2N-2i^)i^2Niy) " V'2^7V-2 (y) V'2]v (2^) 



27r sin^ 6 



sm{hN-2{0 + A9) - f^^m + smihNiO + AO) - (/2iv-2(^^)))] 
N 

[cOs(/2JV-2(e)) COs(/2Ar(0 + A9)) - COs{f2N-2{0 + A9)) C0s(/2Ar (0) )] 



277 sin-^ 6* tan 26* 



N COS 46* 
TT sin^ 6* 



sin A9 



,d f2N{0) 

86 



N sin A9 



2-K sin^ 6* tan 29 



sin A6I 



,g /27v(g) 
961 



(5.29) 



Here, we have neglected the oscillatory term of 0(1) arising from the even function. Finally, we get 



1 



TT sin"^ 29 



(cos46'-cos2 29) 



sin A9 



,d f2NiO) 

89 



1 



Combining all, we get 



^2Ni^^ y) 



= -sin[8iVsin^6iA6i]. 

TT 



3in(a:: ^l^y/{AN ~x)Ax) 
■kAx 



With a; — > 2/, we get the level density 



S'i^{x,x) = i^l^^— ^, 



o<x <m. 



Q<x <AN. 



(5.30) 



(5.31) 



(5.32) 



With Ax 0, (i.e. in the bulk of the spectrum) and r = AxS2^Jj{x, x), we get the "universal" sine-kernel (I5.16P 
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6. UNIVERSALITY FOR SYMPLECTIC ENSEMBLES 

Before we prove universality of the eigenvalue correlation for symplectic ensembles, we would clarify some of the 
confusing notations related to the corresponding SOP. 

SOP corresponding to symplectic ensembles of random matrices can be defined in an interval [xi,a;2] as 



Trf {x)7:'^' (x) - vr^ (x) w{x)dx = / cff i^Wk (^) ^ €' i^W (^) = ^^•>'' (6-1) 



where 



(6.2) 



Here, tt'^'^\x) are SOP defined with respect to w{x). This definition is used in Q and 3 and will also be used in this 
paper to study the statistical properties of the symplectic ensembles. 
An alternative definition is if we write 



(6.3) 



such that 
2 



9j 



7rfix)4'y{x)~n['\x)nfix) w^{x)dx 



X2 



'cj>^f^{x)i,^\x)-<pf{x)i^'\x) 



dx 



4>f {.x)i;'C> {x)dx 



(4), 



(6.4) 



Here ^f^(xi) = 4)f\x2) = 0. The SOP Tif\x) in this definition is defined with respect to uP'ix). This definition is 
used in Q and @ and is used to prove duality between SOP corresponding to orthogonal and symplectic ensembles 
of random matrices. Also, this definition differs from that of [1,0, [2^ for /3 = 4 by a factor 2, which is incorporated 
in the normalization constant. 

We will use (I6.ip and the GCD formula to prove universality. Here we would like to mention that our GCD results 

(4) 

are valid for both these definitions with some minor difference in R (x). 



Jacobi SOP 



We consider SOP (|6.2[) defined with respect to the Jacobi weight (|2.6p . As shown in [1] and 
and T^'j{x) can be written compactly in terms of Jacobi OP P°''^(x): 

T^'2m+l{x) = P2rai^)^ TO = 0,1,..., 

Amix) = P2rl-lix)+mr,iT^'2m-2{x)^ TO=1,2,..., 1t'^{x)^Q, 

where •q2m is a constant, given in Eq. (j6.10p below. On integration, we find the polynomials: 



TTf\x) = T:j{x) 



(6.5) 
(6.6) 



7r2m+i(a;) 

T^2m{x) 



D2m+lP2rn+l{x) + E2rn+lP2rn{x) + F2m+lP2^-i(x) 



{2m + a + b) 
2 

{2m + a + b- 1) 

D2mP2mix) + E2rnP2^-l{x) + -F'2m-P2"m_2 (a^) + ??2m7r2m-2 (a;) , 



m = 0, 1, . 



TO = 0, 1, 



Here, P°' (x) = for negetive j. In (|6.7p and (|6.8p we have used the indefinite integral 

lU + a + b) J Pf{x)dx = P;^,''''-\x), 

= D,+,P;i\ix) + E,+,Pf{x) + F,+,p;i\ix) 



(6.7) 



(6.8) 



(6.9) 
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The integration constants have been put equal to zero because of skew-orthogonahty with tti (a;) . The constants Dj , 

J4) 



Ej, Fj, rjj and g)j are given by 

{i + a + h){i + a + h-\) 



D, 



(4) _ (4) 
92m — 92m+l 



For large j and large m, 



{2j + a + b){2j + a + b-l)' 
ij + a-l)ij + b-l) 
{2j + a + b-l){2j + a + b-2) 



[a-b){3 +a + b-l) 



^ {2j + a + b){2j + a + b-2y 

{j+a-l){j + b-l){2j + a + b-b) 



{j-l){j+a + b-l){2j + a + b-l)-' 



Am + a + — 1 ' 

2a+b+2p(2m + a + l)T{2m + b+l) 
(4m + a + 6 + 1)(4to + a + 6 - l)r(2TO + l)r(2m + a + 5 + 1) ' 



r;, = l + 0(j-^), 



(4) 



{a-b) 

2a+b 



Am? 



+ 0(m-^), 



and in the same approximation, 



Am 



Am 



(6.10) 
(6.11) 

(6.12) 
(6.13) 



^2ra+l{x) = ^[P2m+l{x) - F^^'-ll^)], ^2m{x) = -^[P^™!^) + 2(''+^)/2 (z«,,, (x) ) "^Z^] . (6.14) 



Here in (j6.14p . the non-polynomial term on the right hand side is the large-m approximation for the lower-order terms 
in the series in (|6.8p and has been verified numerically. 

For large to, we use ()6.14|) and the asymptotic formula for Jacobi OP (|4.3p to obtain asymptotic formula for Jacobi 
SOP. For a, b arbitrary and real, e a fixed positive number, we have for 



where 



X = cos 



e < 6I< TT - e, 



(52m) ' J Wa,b{x)TT2,n+l{x) := <t>2l+l(.x) = - 



l[/2,n(^)]+0(TO-3), 



(52m) ' \l Wa,b{^)T^2m{x) -.^ 4>2i{x) 



'1/2 



(4) 



1 



%/ TTTO sin 9 



3[/2m(^)] + l 



0{m- 



4'2+l(^) = 2,/^IcOs[/2m(0)] + 0(to-^: 



4'i(^) 



[sin[/2m(^^)] + l] + 0(TO-2). 



/2m (^) 



2to 



a + b+l 



"+2 2- 



(6.15) 
(6.16) 

(6.17) 
(6.18) 

(6.19) 



Equation (|6.17p is derived from (|6.5p and (|6.18p is obtained by differentiating (|6.16p . Here, we would like to mention 
that to calculate level density and two-point correlation function for the Jacobi symplectic ensemble, (|6.17p and (|6.18p 
are not needed. However, they are important to define the SOP and hence included for completeness. 



B. Level density and 2-point correlation for Jacobi symplectic ensemble 



In this subsection, we calculate the level density of Jacobi symplectic ensemble and show that in the bulk of the 
spectrum, the scaled or "unfolded" |2^ kernel function (S'2^(a;, j/)/S'2^(?/, y)) is stationary and universal. 
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As suggested in Ea. (|2.3ip . to obtain the kernel function, we need to calculate -R2W-1 27V+i' ^'Cn-i 2N^ ^-in- 



?(4) 



?(4) 



2,27V 



and ^^2^-1 2N- I^or this, we use (|2.18p 



2m+2 



j=2m-2 



to get -P2m+i,2m+2- I* is given by 



R 



(4) 



2m+l,2m+2 



A 



92m+2 (2m- + a + b + l)Q2m,2m+lQ2m+l,2m+2 



(4) 
52m 



2D 



2m+2 



where Qj k and Dj are given in the table and ()6.10|) respectively. 
Similarly, we use (|2.19p : 

, 2m+2 , 2m+3 



j=2m-2 



j=2m-2 



to get ,2m+2i -^2m+i,2m+3 ^^^^ -^2m+i,2m+2- -f'or large m, wc get to thc leading order 



(6.20) 



(6.21) 



(6.22) 



R 



(4) 

■2m+l,2m+2 



\ 



52m+2r(2TO + a + 6 + 

l)-E'2m+3Q2m,2m+lQ2m+l,2m+2Q2m+2,2m+3 



(4) 
92 m 



-{a-b) 




{2m 



2D2m+3D2m+2 
b+l) D2rn+1 



2{2m + a + b) L»2m+2 
for a ~ b. 



2m+l,2m+2] + (« - b)0{l/m) 



(6.23) 
(6.24) 



-"-2m+l,2m+3 



52m+2 (2w + a + 6 + 2)(52m,2m+lQ2m+l,2?Ti+2Q2m+2,2m+3 



(4) 
52m 



2D 



2m+3 



o(i), 



(6.25) 



(4) 

•2m,2m+2 



52m+2 (2"1 + a + 6 + l)Q2m-l lQ2m+l,2m+2 



(4) 
52 m 



2£'2m^ 



Oil). 



where Qj^k, Dj and Ej are given in the table and (|6.10p respectively, 
we 

P2N-l,2N-^-N + 0{l), 



For large j, we use Qj,j+i — 1/2. Also for rn = — 1, large N, we have (52tl/52?v-2) — -'^ ^^^d 



R 



(4) 



A^ 



•2iV-l,27V+l 

^-^ + o(i), 



R 



(4) (^-«) I ..^rnv-^ r(4) 

2Ar-l,27V ^ +1" U)(J{I\ ), ^2N-2.2N 



(6.26) 



A^ 

-y + 0(1).(6.27) 
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Finally using (|6.15p and (|6.16p . defined in the 6 interval [e, tt — e], and (|6.27|) in the GCD formula (|2.3ip . we get for 



y — cos 6, X = y + Ay 
{y-x){l-y')si''jl{x,y) 



cos^ 

?(4) 



Ad), 



y — X ~ AO sin 0, 



R 
+R 



■2N-2,2N 
(4) 



nNi^mN-iiy) - nNivmN-ii^) 



2N-1.2N+1 



<t>^2N-2{^)4>2N+l{y) - <t>^2N-2{y)<t>^2N+l{^) 

41v-2(2;)027v(2/) - 41v-2(2/)41v(2;) 



^2N-i:2N ■^^2N-1,2N 



— [cOs(/2Ar(0 + A6)) s\n{hN-2{e)) ~ COs(/2Ar(0)) sin(/2N-2(^^ + A9))] 

47r 

+ ^ [cOs(/2Ar_2(e + AO)) Sm{hNm - COsihN-2iO)) sin(/2Ar(0 + A^))] 



COS 6* 
An sin 6 



[cOs{f2N-2{d + Ad)) COS(/2W(0)) - COs(/2iV-2 (0) ) C0s(/2Ar(^? + A0))] 



[sin[f2N-2{e) - {f2N{e + AO)] + sin[f2N{e) - f2N-2{9 + A6)]] 

47r 



cos6'sin(26') 



47r sin i 

1 

Ait 

cos6'sin{26') 



sin[(2iV 



47r sin 6 



'M^0f2N 

a + b+1 
sin[(2A^ + 



)A9 - 29] + sin[(2iV - 



b+l 



)A9 + 26] 



'-)A9] 



1 

2^ 



sin[2iVA6l] [cos(2e') - cos^ 9] , 



where in the second step, we have dropped 0{N ^) terms. This gives us 

(4) _ sin(2iVAg) 

^2N[^,y) - 2nAes\n9' 

sin(27V(l - y2)-i/2Ay) 



With Ay — > 0, we get the level density 



27rAy 
N 



|2/I<1- 



Ti'Vi ~ y^ 

In the bulk of the spectrum, we get the "universal" sine- kernel 

S^Ni^^y) _ sin27rr 



|2/|<1- 



stN^y^y) 



2nr 



r = Ay5^^(y,y). 



(6.28) 



(6.29) 



(6.30) 



(6.31) 



C. Associated Laguerre SOP 



Now we consider SOP (|6.2p defined with respect to the associated Laguerre weight (|2.7p . As shown in Q and [3], 
Tr'^'^\x) = TTj{x) and TTj{x) can be written compactly in terms of associated Laguerre OP l'^°^\x): 



Am+i{x) = 4m (a;), m = 0,l,. 



7^2m(a;) = ^2^-1(2;) + 



2m + a — 1 



4m-2(a;), m==l,2,. 



2m- 1 

On integration, we find: 

TT2m+l{x) = -i2m+l(2;) +-^2m(a;), m==0, 1, 

71-2m(a;) = -L'2^[x)+L}2^_i[x)+ [— — 7r2m-2, TO = 0,1,. 

\ 2m — 1 



4{x) = 0. 



(6.32) 
(6.33) 

(6.34) 
(6.35) 
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For a = 0, (|6.341 [6735)1 give back the results of [29|, with the observation that any multiple of TT2m{x) can be added 
to T^2m+iix). The normalization constant is given by 



(4) _ (4) 
92m — 92m+l 



' (a) 
'■2m- 



(6.36) 



The results (|6.341 16.35P derive from (|6.32( I6.33P from the indefinite integral, 



the constants of integration in (|6.34p , (|6.35p being zero on skew-orthogonality with tti (x) . 

To obtain asymptotic formula for associated Laguerre SOP, we use (|6.34p . (|6.35p and the asymptotic formula for 
associated Laguerre OP (|4.4p and (|4.5p . To avoid 9 floating inside the argument, for a given x, we choose 9 which 
effectively corresponds to j = 2m + 1/2 in (j4.4p . For 'a' arbitrary and real, e a fixed positive number, we have for 



with 



'a'^^W-''^ (x) 



' 92m''p2m+li^) 



im + 2a + 4) cos^ 9, e < < 7r/2 - ern ^ ^ 



72m+l/2; 



(2m) 



a/2 



1 



2^/ 



(2m) 



a/2 



V mil tan 
= (2to)"/2 



7rm cos sm 

0(1) 



:COs[/2™(0)] + l + 



0(1) 



m^Jx 



sin[/2m(6')] 



/mx 



1 



8y/nm cos-^ 6'sin0 



sin[/2m(0)] 



0(1) 



(2to) 



a/2 



1 /tan0 0(1) 

O V cos /2„i 61) + j= 

2 V mTT mwx 



(6.37) 

(6.38) 
(6.39) 

(6.40) 



y 92m,±2V2m±2\^ ) " ^ 



1 



2 V TTTO COS sin"^ 9 



COs[f2rn{9) T 49] + 1 + 



0(1) 



Similar relations hold for (jOS)) . (jOOl) and fOO)) as 



± 1. Here 



/2m(6') = (2m + 1 + a/2)(sin20 - 29) + 



(6.41) 



(6.42) 



In deriving (|6.37p we have used the large-m approximation 



2to 



7r2r 



= - {L^2l{^) + ~ I (2m)'^/^(^.(x))-/^ 



which follows from the three-term recursion and a sum rule for L, (x) [1 



(6.43) 



D. Level-density and 2-point correlation for associated Laguerre symplectic ensemble 



In this subsection, we calculate the level density of associated Laguerre symplectic ensemble and show that in the 
bulk of the spectrum, the scaled or "unfolded" [1^ kernel function (S'2^(x, 2;)/52'^(i/, y)) is stationary and universal. 

To study the kernel function (|2.3ip . we need to calculate R^^_i 2n+1' ^2JV-i 2^' -^2lv-2 2N ^^"^ ^2N-i 2N- ^^^r 
this, we use (|2.18p and (|2.19p for associated Laguerre weight: 



2m+2 



j=2m-2 
2m+3 



^ (^2.'^2m-|-l(^) — E/ ^2m+l,j't''j \^ 
j=2m-2 



2m+2 



/ , -''-2m 
j=2)ri-2 



(6.44) 
(6.45) 
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From which, using properties of OP, we get 



d(4) 

2m+l,2m+2 



\ 



(4) 
g2m+2 

(4) 
92m 



(m+1), R2m+l,2m+3 



(4) 

32771 + 2 

(4) 
92m 



(to+ l)(2m + 3). 



(6.46) 



R 



(4) 

■2771,2771 + 2 



A 



(4) 
ff2m+2 

(4) 
92m 



m+l)(2m + l), 4t',+i,2,„+2 



\ 



(4) 

ff2771 + 2 

(4) 
92m 



(to + l)(4m + a + 4). (6.47) 



For m = — 1, large TV, we have {,g2N / 9^n -2) — 1- Finally for large N , the recursion coefficients are given by 



p: 



(4) 



27V-l,2Af 



TV + 0(1), 



R 



(4) 

2Af-l,2Af+l 



-2iV2 + 0(iV), 



(4) 

2N-2,2N 



2N^ + OiN), Ri%_,^^j^^m^ + 0{N). (6.48) 



Also, since 9^il^2/ 9tn — Ij we can write 



(4) (4) 
92N-292N 



w{x)w{y)TT2N{x)TT2N-l{y) 



(4) 
5'2Ar 



V w{x)w{y)'K2N{x)'n:2N-iiy) 



■ 7J— Vw(a::)w(2/)7r2Ar (a;)7r2Ar_i (y) 

'''2N 



-i2N) "\Jw{x)'w{y)'K2N{x)TT2N~i{v)- 



(6.49) 



Finally using the asymptotic results (|6.38p . (|6.37p . (|6.4ip . defined in the 9 interval [e,7r/2 — eN 2]^ and the large 
results (|08l) and (fOO]) in the CCD formula ([23T|) . we get for 



y = (87V + 2a + 4) cos^ 



y + Ay, 



y(y- a;)5y(a;,y) = R2n-2,2N 4w (2^)4^-1 4^(^)4^-1 W 



-i?: 



^''^ 4w-2(2;)4w+i(2/) - 4w-2(2/)4w+i(-^) 

^P2N~1,2n) 41v-2(2^)47v(2/) - 47V-2(2/)4lv(2^) 
{sin[/2Ar_2(^^) - f2N{d + M)] + sin[/2Ar(0) - f2N^2{0 + AO)]} 



2N-1,2N+1 
\^2N-1,2N 

N 

27r sin^ 

TV cos 20 {cOs[f2N-2{0 + AO)] COs[(/2Ar(0)] - COs[/2^-2(e)] COs[/2Ar(0 - 

47r COS sm 



NCOS 49 iV sin 46* cos 26* 



TT sin 6 A-K cos sin"^ 9 



sm I „ „ AO 



= SNcos^e 



[cos46l-cos2 261] 
27r sin^ 29 



sm 



961 

df2N 

do 



AO 



(6.50) 



which gives 



Thus we get 



^2Ni^' y) 



sin( iy-VsiV-yAy 
27rAy 



< y < 87V. 



(6.51) 
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With Ay 0, we get the level density 

Si%,y) = 0<2;<8iV. (6.52) 

Finally, in the bulk of the spectrum, we get the "universal" sine-kernel (|6.3ip . 

7. CONCLUSION 

In [| and the authors prove "universality" for the entire Jacobi ensembles of random matrices using SOP, 
which were written in terms of OP. The asymptotic properties of OP were used to obtain asymptotic formulae for 
SOP. Finally, the summation in {x, y) was replaced by an integral for large iV. However, using asymptotic results 
of the SOP to calculate the kernel function, before completing the sum lacks mathematical rigor. 

In this paper, we have shown that SOP <i>lf''(a;) and '^\i\x) corresponding to classical weight satisfy three-term 

recursion relations in the 2x2 quaternion space. Using this, we obtain the kernel functions s'^^{x^y), /3 = 1,4, 
for the entire family of finite dimensional Jacobi ensembles of random 27V x 2N matrices. As iV ^ cxd, we use the 
asymptotic results of the SOP in the range [xi + e,X2 — over which they are defined, to prove that in the bulk of 
the spectrum, the correlation functions are universal. 

Here, we would like to mention that our GCD results are valid in the entire complex plane. Hence to study 
statistical properties of the eigenvalues of orthogonal and symplectic ensembles away from the bulk, we need to use 
Plancherel-Rotach type formula for these polynomials defined outside the range [xi -I- e, X2 — e], something which has 
already been done for the unitary ensembles. We wish to come back to this in a later publication. 

We would like to emphasize that the key step in deriving the asymptotic results of SOP is to obtain and solve finite 
term recursion relations between SOP and OP. Recently, this technique has been used to obtain bulk asymptotics 
of SOP corresponding to quartic double well potential [2]. Till date, this seems the easier method to study the 
asymptotic behavior of SOP rather than solving the 2d x 2d Riemann Hilbert problem (Gl.fisj. 

Finally, a word on duality. Results in section 3 makes us wonder if the two families of SOP are really different or 
there exists a simple mapping between them. We have seen the existence of a simple relation (|3.14p between SOP 
arising in the study of the associated Laguerre ensembles. The other ensembles do show similar pattern, although we 
are unable to come out with a general formula. A deeper theoretical understanding is needed to obtain the mapping 
between these two families of SOP. 
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